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ABSTRACT 

We investigate how inhomogeneous quintessence models may have a specific signa- 
ture even in the hnear regime of large scale structure formation. The dynamics of 
the collapse of a dark matter halo is governed by the value or the dynamical evolu- 
tion of the dark energy equation of state, the energy density's initial conditions and 
its homogeneity nature in the highly non-linear regime. These have a direct impact 
on the redshift of collapse, altering in consequence the linearly extrapolated density 
threshold above which structures will end up collapsing. We compute this quantity for 
minimally coupled and coupled quintessence models, examining two extreme scenar- 
ios: first, when the quintessence field does not exhibit fluctuations on cluster scales and 
below - homogeneous dark energy; and second, when the field inside the overdensity 
collapses along with the dark matter - inhomogeneous dark energy. One shows that 
inhomogeneous dark energy models present distinct features which may be used to 
confront them with observational data, for instance galaxy number counting. Fitting 
formulae for the linearly extrapolated density threshold above which structures will 
end up collapsing are provided for models of dark energy with constant equation of 
state. 
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1 INTRODUCTION 

Measurements of the luminosity-redshift relati onship from 
observations of supernovae of type la (SNIa) iRiess et alJ 
1200 ll: IPerlmutter et alj Il999l : iKnop et alj|2003l) . the mat- 
ter power spectrum of large scale structure as inferred from 
galaxy redsh ift surveys like the Sloan Digital Sky Survey 
(SDSS) (iTegmark et alj|2004) an d the 2dF Galaxy Redshift 
Survey f2dFGRS) JColless et alJll9^8^ ■ and the anisotropies 
in the Cosmic Micr owave Background Radiation (CMBR) 
l|Spergel et alJl2003^ . reveal that a mysterious constituent 
with negative pressure, so-called dark energy, accounts for 
seventy percent of today's mass-energy budget and is caus- 
ing the expansion of the universe to accelerate. 

Despite of its major importance in explaining the astro- 
physical data, the nature of dark energy is one of the greatest 
mysteries of modern cosmology. The simplest candidates for 
this unknown ent i ty inc lude the cosmological constant (see 
e.g. ICarroll et al.| (|l992h ') and a scalar field, commonly des- 
ignated quintessence, which varies across space and changes 
with time (see e.g. iRatra fc Peebles (■198&) '). 
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Though the nature of dark energy is unknown, one can 
still try to infer its properties from its effects on cosmic 
structure formation. In fact, the behaviour of linear per- 
turbations in a scalar field and its effect on large struc- 
ture formation has been investigated by a number of au- 
thors (see e.g. (iFerreira fc Jovcgll998^ '). However, the be- 
haviour of quintessence during the non-linear gravitational 
collapse is not well understood and is currently under in- 
vestigation (see e. g. Wottcrich (2002); Amcndola (200;_ 
iMota fc van de Bruck ( 20041) : iMaor fc Lahavl ll2004l) : IWam 
|2005j); IPercival (i2005i) ). Usually, it is assumed that the 
quintessence field does not exhibit density fluctuations on 
cluster scales and below. The reason for this assumption is 
that, according to linear perturbation theory, the mass of 
the field is very small (the associated wavelength of the par- 
ticle is of the order of the Hubble radius) and, hence, it does 
not feel matter overdensities of the size of tenth of a Mpc or 
smaller dWang fc SteinharddHool) . 

The assumption of neglecting the effects of matter per- 
turbations on the evolution of dark energy at small scales 
is indeed a good approximation when perturbations in 
the metric are very small. However, care must be taken 
when extrapolating the small-scale linear-regime results 
to the highly non-linear regime. Then, locally the flat 
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FRW metric is no longer a good approximation to de- 
scribe the geometry of overdense regions. Highly non-linear 
matter perturbations could, in principle, modify the evo- 
lution of perturbations in dark energy considerably, and 
these could, in turn, back-react and affect the evolution of 
matter overdensities. Moreover, it is natural to think that 
once a dark matter overdensity decouples from the back- 
ground expansion and collapses, the field inside the clus- 
ter feels the gravitational potential inside the overdensity 
and its evolution will be different from the background evo- 
lution. This is indeed a general feature of not only cos- 
mological scalar fields whose properties d epend on the lo- 
cal density of the region they " li ve in" iMot a fc Barrow! 



2004al: Ixhourv fc Weltmanl 120041: iMota fc Barrow 20041: 



Khourv fc Weltmanl l200A- Iciifton. Mota fc Barrow! 120041) . 

but also of massless particles, such as, photons (example 
of astrophysical effects are: Lensing, Riess-Schiama effect, 
Sac hs- Wolfe effect). 

iBean fc Maeueiiol i2002h suggested that the 
quintessence field could have an important im p act i n 
the highly non— linear regime. IWetterichI ll200ll . l2002l) : 
lArbev et all (j gOQl) noted that the quintessence field could 
indeed b e important on galactic scales . It was put for- 
ward by iGuzman fc Urena-Lope j (|2003); lAlcubierre et alJ 
i2002h that it could in fact be responsible for the ob- 
served flat rotation curves in ga l axies. Other authors 
JPadmanabhan fc Choudhurvl 120021: IPadmanabhaiil 120021: 
iBada et al.l 120031: ICaussd ' l200^r (iiscussed more exotic 
models, based on tachyon fields, and argued that the 
equation of state is scale-dependent. 

If it turns out that the effects of dark matter den- 
sity perturbations and metric influence perturbations of 
quintessence on small scales, this could significantly change 
our understanding of structure formation on galactic and 
clus ter scales. 

iMota fc van de Bruckl (^2004) have shown that proper- 
ties of halos, such as the density contrast and the virial 
radius, depend critically on the form of the potential, the 
initial conditions of the field, the time evolution of its equa- 
tion of state and on the behaviour of quintessence in highly 
non-linear regions. In reality, the dependence on the inho- 
mogeneity of dark energy is only important for some dark 
energy candidates. If the dark energy equation of state w is 
constant, the differences between the homogeneous and in- 
homogeneous cases are small, as long as the equation of state 
does not differ largely from w = —1 iMota fc van de Bruckl 
l2004l:lMaor fc Lahavll2o'oi) . 

If dark energy is indeed a non-negligible component in- 
side matter overdensities, then the dynamics of structure 
formation may be strongly dependent on the nature of dark 
energy. Different dark energy models will contribute in very 
particular ways to the gravitational potential of the over- 
density via the Poisson equation. This may result in specific 
signatures in the formation of large scale structures. 

In this paper, we investigate how inhomogeneous 
quintessence models may have a specific signature even in 
the linear regime of large scale structure formation. In par- 
ticular, we investigate how the time of collapse is affected by 
the inhomogeneity of dark energy and how this is refiected 
on the linearly extrapolated density threshold above which 
structures will end up collapsing, i.e. 5c{z) = 5l{z = Zcoi)- 
This work extends upon previous studies in that we examine 



the evolution of matter overdensities as a function of a time- 
varying dark energy equation of state and its homogeneity 
nature in the non-linear regime. 

The article is organised as follows. In section 2 we inves- 
tigate minimally coupled dark energy models and describe 
briefly the spherical collapse model and its dependence on 
the homogeneity nature of dark energy. In section 3 we gen- 
eralise our analysis to the case of coupled quintessence mod- 
els. A summary of our principal results is given in section 



2 MINIMALLY COUPLED DARK ENERGY 
MODELS 

We will consider a spatially flat Friedmann-Robertson- 
Walker Universe with scale factor a{t). The cosmic dynam- 
ics is determined by a background pressureless fluid (cor- 
responding to dark and visible matter), radiation and dark 
energy. The governing equations of motion are 



H = 



(1) 

(2) 
(3) 

(4) 



+ Pb+ P4> +P<t>) . 

p^ = -3H{p^+p^), 
subject to the Friedmann constraint 

■ff^ = Y (pfl + P4,) ■ 

where H is the ratio of expansion of the Universe H = d/a, 
and — 8tvG. Pb and Pb are, respectively, the energy den- 
sity and pressure of the background fluid (dust and radi- 
ation). In this work we consider two possible scenarios for 
the nature of dark energy. If the dark energy is a perfect 
fluid then its energy density and pressure are related by the 
equation of state p^ = w^p^ and p^, = Q.^opo/a^'^^'''^^^ ■ Al- 
ternatively, dark energy can be described by a dynamical 
evolving scalar field rolling down its potential V{(j)). In this 
case, it is defined energy density and pressure of a scalar field 
as, p^ = (fy^ /2 + V{(j)) and = (/>^/2 — V{(l)), respectively. 
The equation of motion for the scalar field is, 

if} = -3//0 - ^ . (5) 



2.1 Dark energy models 

In this work we will be exploring dark energy models for 
which = —1 (the cosmolog ical const a nt), = —0.8, 
ui0 — —1.2 (phantom energy, ICaldwell ll2002l) 'i and two 
cases where the dark energy is the result of a slowly evolv- 
ing scala r field in a po t ential with two exponential terms 
(2EXP) jBarreiro et alJ J2000h ) 



V{4>) = Vo (e 



(6) 



We have chosen the pair (q,/3) = (6.2,0.1) (a) and (a,/3) = 
(20.1,0.5) (b) as they both provide an equation of state at 
present w^o — —0.95 though having distinct evolutions at 
higher redshift. The equation of state for (a) approaches zero 
at a higher rate than for (b). What mainly distinguishes 
these two models is the contribution of the field at high 
redshift. As can be seen from Table. (a) provides a contri- 
bution of dark energy that is non negligible at high redshifts. 
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Table 1. Evolution of the effective equation of state and con- 
tribution of the dark energy energy density with rcdshift for two 
pairs of parameters of the 2EXP model. 



We expect, therefore, that this model will provide features 
with a stronger departure from a pure cosmological constant 
than the remaining models. 

One should also emphasize that the values of the present 
value of the equation of state and its running (wq and 
w'o — dw/dz{z = 0)) of these models are within the current 
bounds dete r mined from supernovae la (SNIa) observations 
jRiess et al.l ^00^). Models (a), w = -0.8 and w = -1.2 
correspond to limiting cases consistent with these observa- 
tions. We have chosen them to estimate the largest range of 
departur es on 5c from the cosmological constant case. See 
however, Ijassal et al.l ^2005^ where it seems that model (a) 
is not consistent with WMAP data. 

In Table Q we show how the effective equation of state 
(i.e. weighted average of the equation of state between red- 
shift and z) depends on redshift for these models. We are 
assuming that for all of these models n^o = 0.7. 

Other mode ls of dark energ y such as the in- 
verse power law dZlatev et alJ (Il999h ') or SUGRA models 
iBrax fc Martini Jl999l) ) have negligible contribution of the 
scalar field at high redshift, hence, they are not entirely dis- 
tinct from a Wcj, — constant dark energy and we will not 
consider them here. 



2.2 Spherical collapse model 

In the remaining of this section we will be using the spherical 
collapse model to describe the gravitational collapse of an 
overdense region in minimally coupled dark energy models. 
The radius of the overdense region r and density contrast 5 
are related in this case by 1 -I- 5 = Pmc/pm = [a/r)^, where 
pmc and pm are the energy densities of pressureless matter 
in the cluster and in the background, respectively. In the 
next section we will study the spherical collapse model for 
coupled quintessence models in which case this relation is 
modified. 

The equation of motion for the non-linear evolution of 
the density contrast is 



+ — [(1 + 3w^Jp^^ - (1 + 3w^)p^] (1 + 5) , (7) 

where we have considered the possibility that dark energy 
also clusters, i.e. p^c 7^ p^. In general, the evolution of Pde i n 
the cluster can be written as jMota fc van de Brudd (12004^ ') 



(8) 



where describes the dark energy loss of energy inside the 
cluster. Note that r/r is given by 



raid 
r ^ a ~ 3 1+5 



(9) 



hence, the syst em of equations closes. 

Following jMota fc van de Bruckll200il . we study the 
two extreme limits for the evolution of dark energy in the 
overdensity region. In the first we assume that dark energy 
is homogeneous, i.e. the value of p^ inside the cluster is the 
same as in the background, with 



+ 1 



J: 



(10) 



in which case Eqs. ^ and ||SJ| are equivalent. In the second 
limit, dark energy is inhomogeneous, collapses with dark 
matter such that F^ = 0. In this case, the equation of motion 
for the scalar field inside the cluster is 



= — - 
r 



dF(0c 



(11) 



Another form of F ^,, with a different limit and situat ion, was 
also considered in (|Maor fc Lahavl20o3 : IWangl200,'J l. but we 
will not consider it here. 

The linear regime of Eq. (|7J defines the linear density 
contrast 5l determined by the equation 



5l = —2H5l + — lp,n5L + (1 + 3w^)p^ 5^ + 3p4,5w^] 



(12) 



where we have defined 5^ = 5p^/p^. In the homogeneous 
case, 5^ — 5w^ — 0. On the other, for inhomogeneous dark 
energy with a constant equation of state, 5w^ = and 5^ 
satisfies the equation of motion 



64, = -~2H5^ + -2"(1 + [pm5L + (1 + Sw4)p454,] 

When dark energy is sourced by a scalar field, we have 

dV . , 



(13) 



(14) 



1 dV r 



where the perturbation of the field, 5(j>, satisfies 
54> = -SH5^ - ^ 5(j, + 5l <i> . 



(15) 



(16) 



Equations I12II and 1151 1 are equivalent to the ones found in 
iHwang fc Noh ( 2001). 

In the case where the equation of state of dark energy is 
a constant, one can derive that when F^ = the non-linear 
evolution of the energy density of dark energy inside the 
collapsing region relates to the the one in the background 
through 



1 + ^ 

1 + 5,. 



P0 : 



(17) 



for some initial perturbation 5i. In Fig.Qwe compare the de- 
pendence of the density contrast and linear density contrast 
with redshift for a homogeneous and inhomogeneous scalar 
field dark energy model. We see that in this case the density 
contrast and the linear density contrast evolve with different 
slopes from the homogeneous to the inhomogeneous cases. 
This is a result of a non negligible contribution of the scalar 
field at high redshifts in this particular model. In the case of 
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Figure 1. Evolution of tlie density contrast 5, and linear den- 
sity contrast 5^ for the 2EXP model with (a,/?) = (6.2,0.1). 
Solid line: density contrast in homogeneous case; dashed line: lin- 
ear density contrast in homogeneous case; thin solid line: density 
contrast in inhomogeneous case; thin dashed line: linear density 
contrast in inhomogeneous case. 



a dark energy model of constant equation of state, the var- 
ious evolutions are essentially indistinguishable except for 
the later stages when the dark energy starts to dominate at 
low redshift. In general terms, differences between homoge- 
neous and inhomogeneous models depend on the equation of 
state of dark energy (its value and if it is dynamical or not) 
and on the contribution of dark energy for the total energy 
budget of the Universe at high redshifts. 

We can now compute the linearly extrapolated den- 
sity threshold above which structures will end up collaps- 
ing, i.e. 5c{z) — 5l{z = 2coi). This quantity, is essential 
to compute the number of col lapsed structures followin g 
the Press-Schechter formalism jPress fc Schechteil il974l) '). 
In our analysis we integrate the equations of motion from 
very high redshifts {1 + Zi = 10*) to minimize the effect of 
the decaying mode on the evolution of the perturbations. 

For homogeneous scalar field dark energy, the solution 
of Eq. 1121 for a component (/> of fixed contribution to the 
total energy density fi^ = 1 — Sim = (as it is the case 

of dark energy models with an exponential term dominant 
at early times) and neglecting the radiation component, is: 



4(..o4(i+A-)(^) 



(A-l)/4 



(18) 



where A = ^/25 - and we have neglected the de- 

caying mode. Given that 72/q:^ — 24(1 — Qm) and that 
1 + z (X t~^^^, we obtain the time dependence S oc t"*. where 
m = (Vl + 24r2m-l)/6, found bv lBagla et"ai] lj200,'j) . how- 
ever, in the context of a scaling tachyon^. In an Einstein-de 
Sitter Universe (i.e. a — oo), solving Eq. with 1+Zi = 10* 
for an initial perturbation Si — 10~^ , the overdensity re- 
gion coll apses at redshif t z^oi — 2.558, hence Sc = 1.686 
(see e.g. IPadmanabhan I (Il995ll ). This classical value was 
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Figure 2. Linear density contrast at the time of collapse in 
the homogeneous (top panel) and inhomogeneous (bottom panel) 
dark energy scenarios. 



used in the seminal paper of iPress &: SchechteJ (^7^. For 
a model with one pure exponential potential with a = 6.2 we 
would obtain for the same initial conditions ^coi = 0.522 and 
5c = 1.681. This is the asymptotic value of Sc for the dark 
energy model considered, (a, /3) = (6.2,0.1), at high red- 
shifts (see upper panel of figure |5J. At low redshifts, when 
the dark energy component becomes important, the value of 
Sc decreases with the redshift of collapse. 

The computation of the value of Sc at the time of col- 
lapse for inhomogeneous dark energy must consider in addi- 
tion the terms in S^ and Sw^ in Eq. (1121 . We can set them 
to vanish as initial conditions, however, they are going to 
evolve in cosmological times, in fact, for a pure exponential 
potential, Sl and 5^ will approach an attractor solution. In- 
deed, defining the dimensionless quantities: 

_ k5(I>' S'r 



IT 



-'L 

Sl 



(19) 



where a prime means differentiation with respect to In a, we 
can rewrite Eas. 1121 and 1161 in the form of a system of first 
order differential equations 



9 , 3 3 

9 6 1 

Za ot 2 



2^3 
T + 



3/ 3 
2 V a 



This system has stable critical points at 

_3_ _ _3_ 

7a ' 7a 



tpc = 



Tc^l.. 



(20) 
(21) 

(22) 
(23) 



where we have used the well known scaling solution result 
for an exponential potential Kcj)' = 1/a and k^V — 9H^ /2a^ 
jCopeland. Liddle fc Wandslll99^ . Substituting back into 
S^ and Sw^ we obtain that S^ = —Sl/U and Sw^ = 5(5l/14 
independent of a. Upon substitution into Eq. 1121 we find 
the interesting result 



Sl = -2HSl 



Ih^Sl.. 



(24) 
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i.e., the linear overdensity evolves in the attractor as if there 
is only dust in the Universe. However we should not expect 
to obtain 5c = 1.686 as in the Einstein-de Sitter Universe, 
for two reasons. Firstly, the system takes a finite amount 
of time before reaching the attractor solution. Secondly, the 
equation of state evolves inside the collapsing region from 
ui = at high redshifts to ui = 1 at the time of collapse 
(see Fig. thus increasing in the cluster favoring the 
clustering of dark matter. Hence, for the same initial condi- 
tions, the collapse occurs earlier, when the growth factor is 
smaller and therefore 5c < 1.686. 

One can understand why the equation of state of dark 
energy evolves inside the overdensity region by inspecting 
equation (IIH . Indeed, after the turn around r/r becomes 
negative switching the frictional into an anti-frictional term 
in the equation of motion of the scalar field. As the scalar 
potential becomes less and less important the kinetic energy 
of the field approaches the asymptotic evolution (j>^ oc 
(therefore an effective w^^ = 1) which would in principle 
overtake the energy density of pressureless matter as pYnc oc 
. In reality, this will never happen, since the singularity 
in the collapse will not be reached. The dark matter halo 
will virialise much before and a dynamical equilibrium will 
be reached, where the halo has a constant radius (the virial 
radius). 

It is important to notice that though the equation of 
state for dark energy inside the overdensity can become pos- 
itive, it is still negative in the cosmological background (see 
figure]^. Hence, the Universe's background expansion will 
not be affected by the local behaviour of the dark energy 
inside clusters. The usual late time accelerated expansion 
still occurs as normally measured by SNIa. The only effects 
due to the positivity of the dark en ergy equation of state 
occur only inside the overdensities iMota fc van de Bruckl 
I2OO4I) . Moreover, within the models investigated, the dark 
energy contribution inside the cluster is always subdomi- 
nant when compared to the one of dark matter and it only 
becomes dominant very near the collapse. We should not 
expect, therefore, any unusual effects on the virialisation 
process resulting from the positivity of the equation of state 
inside the cluster, at least for those models with negligible 
background dark energy contribution at hig h redshifts. Sig- 
natures in lensing and X-ray obse r vations llLopes fc MilleJ 
I2OO4I: iLopes. Mota fc MiUeil I200I ISereno fc Longd |2004^ 
may be noticed though, if dark energy provides a non neg- 
ligible background contribution at high redshifts. 

In Fig. 121 we show the dependence of 5c on the redshift 
of collapse for homogeneous and inhomogeneous dark en- 
ergy. Note that for the cases of inhomogeneous scalar field 
dark energy, the value of 5c at high redshifts is consider- 
ably lower than for w^t, = constant cases. It is worth noting 
that for an inhomogeneous dark energy, the evolution of 5c 
with redshift is significantly altered for the 'phantom en- 
ergy', with 5c decreasing from 1.753 to 1.686 for increasing 
redshift. These variations of the density threshold will obvi- 
ously modify the predicted numbers of collapsed objects. 

The reason why one finds a difference in 5c when com- 
paring homogeneous and inhomogeneous dark energy mod- 
els, becomes clear when one analyses the Poisson equation 
for the gravitational potential, $, inside the overdensity. 




Figure 3. Upper panel: evolution of the dark energy equa- 
tion of state inside the overdense region (dashed line) and in 
the background (solid line) for a cluster that collapses today 
(^col = 0)- This figure results from numerically integrating the 
equations of motion of the scalar field with the 2EXP model with 
(a, 0) = (6.2, 0.1). The turn around occurs at 2ta = 0.7 which is 
shortly before the redshift at which the equation of state inside the 
overdensity itself turns its evolution towards unity in opposition 
to its counterpart in the background that continues to decrease 
to more negative values. Lower panel: evolution of the dark en- 
ergy (dashed line) and dark matter (solid line) energy densities in 
the overdense region. The dark energy contribution only becomes 
dominant very close to the collapse when vja,^ approaches unity. 



In the case of homogeneous models 
^2 

= Y [Pm, + (1 + 3™0)/9^J 



= -^\{^+ S)pm + (1 + 3w^)p^ 



(26) 



as pif,^ — p0. However for inhomogeneous models, making 
use of equation 1171 for a constant equation of state, we 
have 



K 

T 



(1 + 5)pm + (1 + 3™^) ( j P4> 



(27) 



= '^{ptotal + 3ptotal) . 



(25) 



Comparing equations (I26II and 12711 one notices that in the 
inhomogeneous case the contribution of the scalar field en- 
ergy density becomes increasingly important as the density 
contrast grows. Since is generally negative, the source 
term in equation 1251 is smaller than in the homogeneous 
case. Conversely, if is more negative than —1 the contri- 
bution of the scalar field is less and less important and the 
source term is larger than in the homogeneous case. In con- 
clusion, inhomogeneous dark energy models will contribute 
to the gravitational potential inside the overdensities in a 
different way to homogeneous ones. These will clearly affect 
the whole dynamics of the formation of structures. In par- 
ticular the late stages of the non-linear collapse of overden- 
sities, such cis th e time of collapse, virial rad ius and matter 
density contrast iMota fc van de Bruckl2004^ . In the case of 
dark energy candidates with a dynamical equation of state, 
such as scalar fields, the effects will be even more interest- 
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ing. In these cases the source term in equation 12511 wiU 
vary in time since w,^^ will change during the formation of 
the collapsing overdensity, resulting in a more complex evo- 
lution for the gravitational potential. Recall that w^^ may 
even become positive and approach w^^ — 1 during the late 
phase of collapse. Hence the dark energy pressure in equa- 
tion (I25II will source the gravitational potential positively, 
in opposition to the usual homogeneous models where 
is always negative. The effects will be stronger at low red- 
shifts when dark energy dominates the universe density in 
the background universe. 

For a constant equation of state w the following expres- 
sions provide an accurate fit to the evolution of the density 
contrast as shown in Fig. |5| 

5,{z) = ^(12^)'/' [l+ax + Px"] , (28) 

where x = log(f2m(z)) and for homogeneous dark energy we 
have 



O.OOeiw' + 0.0327W + 0.0403 , 

2 



(29) 

P = -0.0163w^ - 0.0294™ -0.0118, (30) 

whereas for inhomogeneous dark energy it reads 

a = 0.1198w^ -)- 0.6226w-h 0.5170, (31) 
P = 0.2022w^ + 0.28771U + 0.0860 . (32) 

These expressions are valid for —1.2 < w < —0.6. 
3 COUPLED QUINTESSENCE 

Let us now look at the case when the scalar field has a 
coupling w ith all or part of the da rk matter (see for e.g . 
Amendolal Eo 03): Tocchini-Valcntini fc Amendola I J2002ri : 
Amendola fc Tocchini-Valentini. (.20031 ). In these models. 



inhomogeneities in the quintessence field may appear due 
to two main reasons: The first, which is the same as in min- 
imally coupled models, is the change of the local geome- 
try of the region where the overdensity "lives in", which 
is just the general relativistic effect of the spacetime defor- 
mation. The second is the dragging of the scalar field by 
the dark matter particles caused by the coupling between 
them. This second effect is similar to the one which occurs 
in scalar-te nsor theories and leads to inhomogeneities in the 
scalar field JClifton. Mota fc Barr'owll2004l: iMota fc Barrow! 
Hooii). Due to these causes, it is then natural to expect 
stronger inhomogeneity effects in this kind of models than 
in the minimally coupled one. 

The background quantities and the ones who live inside 
the collapsing region are given by 



PcDM ~ PO^ 



Pun 



PcDMc 



Pof^umO ( — ) I — 
\ai / \ a 

V fli / \ a 

(1 + <5i)P0fiumc0 ^-^^ (^-p 
(1 + Si)pQQ,cDMcO (~ 



i3(«-B(0o) 



(33) 
(34) 
(35) 



(36) 



where the subscripts "um" and "cDM" mean uncou- 



pled matter and coupled dark matter, respectively. Un- 
coupled matter corresponds to both baryons and to 
uncoupled dark matter. The function B{(j)) represents 
the coupling between dar k energy and d a rk m atter. 
In the model di scussed by iHolden fc Wand3 ll200Cfl and 
lAmendolal (l2000t) . B{<P) = -CKfj), where C IS a con- 
stant. But othe r forms for this function ha v e als o 
been suggested tA mcndola fc Tocchini-Valenti'nil i2002l) : 
iMainini fc Bonometo (2004)). 

The total energy densities inside the cluster and the 
background are therefore, pm = Pum + PcDM and pc = Pumc + 
PcDMc which evolve accordingly to 



„a dB 
= —o—pm H — rrPcDM<? 
a am 



Pc 



= — 3-pc 
r 



dB 



-pcDMc<Pc ■ 



(37) 
(38) 



If we take for initial conditions (fic{a, = ai) — (j>{a = Ui) 
we can define ri = Oi/(l -|- 5i)^^^ and then we can write 
1 + 5 = Pc/pm, as 



1 + 5: 



icDMcOe 



1(<I,^)-B(M 



(39) 



It will become useful to write the ratios G — G{(j)) = 
PcDM/pm and Gc = G{4>c) = PcDMc/pm(l + 5) iu the fol- 
lowing form: 



(40) 



ncDMOeSW)-S(^o) ^umO ' 

The time derivative of the density contrast will now 
have a component coming from the coupling dB/d(j> in the 
equations above 

<5 = 3(l + <5) [^-^1 +(l+5)f(0), (41) 
(compare with equation @) where have defined F[(j>) as 

dB PcDMc ; dB PcDM 



{l + 5)F{d)) 



(1 + 5) 



dB 



Gc 



Pin 
dB 



{1+5). 



G<f) 



(42) 
(43) 



and it results that F{<j)) is 



F(0) 



G 



G 



(rB_ 

d^' 

d^B 
d^' 



+ 



+ 



dB 



'<^^(1-G.) + ^, 



dB\ -2 r^^.dB- 
-d^ ^^^-^^+d^' 



(44) 

The equations of motion for the evolution of the scalar field 
inside and the background are in this case: 



^d- dV dB 
= -6-(p - -TT - -JtPcDM 
a dm dm 



T • 

r 



dV 

d(f)c 



dB 

-r— PcDMc H — — 



(45) 
(46) 



Using these equations we are now able to obtain the 
modified equation for the non-linear evolution of the density 
contrast 



5 = 



4 5^ 



-2^ [5-[l+5)F] + '^Pin{l + 5)5+\- 
a '- -'2 3 1 
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Figure 4. Linear density contrast at tlie time of collapse in Figure 5. Ratio S'JSc where prime is difTerentiation with respect 

the homogeneous (top panel) and inhomogeneous (bottom panel) redshift for homogeneous (top panel) and inhomogeneous (bot- 

coupled quintessence scenarios. t""^ P^'^cl) ^^^P'^d quintessence scenarios. 



+ — [(1 + Sm^Jp,^, - (1 + 3w4,)p^] (1 + <5) 
-^FS+l{l + 5)F^ + {l+S)F. 



(47) 



From linearizing the expression of the density contrast one 
gets 

5l = -2H{5l -.f) + f 

+ [PmSh + (1 + 3w</,^)(50 P4, + Zp^Sw^] . (48) 

where 5cj> is determined through the equation of motion 
dB 



■G p^5l + {5l- f)(t) 



and / is the linearization of F: 

2 



/ = G 



dB 



dB 



cfB 
d^ 



(49) 



(50) 



Assuming a pure exponential potential V{(j)) = 
Voexp{aK<j)) with a = 10, and B{(f)) = —Cucj), we com- 
pare in figure 0] the evolution of the linear density contrast 
at the time of collapse, 5c, between homogeneous and inho- 
mogeneous dark energy. We take into consideration the case 
when all the dark matter is coupled with the scalar field, 
or in other words, that the uncoupled matter is only the 
baryons, hence f2umO = ^ho\ and the possibility that only 
a small fraction of the dark matter feels the field in which 
case we have taken fiumo = 0.25 in figure |1| Similarly to 
what happens in minimally coupled models, there are clear 
differences in the evolution of 5c between inhomogeneous 
and homogeneous models. As expected, these differences are 
more distinct in the case where the amount of coupled dark 
matter is larger (dashed-lines in figure 2J- 

A particular feature is the oscillating behaviour of 5c 
that is amplified in the inhomogeneous case. These oscilla- 



tions are a result of the oscillations performed by the scalar 
field around the late time attractor solution. As we can see 
from Eqs. 1341 and 1361 . these oscillations induce a corre- 
sponding oscillation in the PcDM and pcDMc components (see 
Fig. 7 of Copeland c t al. (20Qi))- 

Though these oscillations appear inexistent in the ho- 
mogeneous scenario they are nonetheless present as can be 
seen by computing 5'c/5c = {d5c/dz)/5c, see Fig. El They 
only seem inexistent because the mean 5c evolves quicker 
than the amplitude of the oscillations. For the inhomoge- 
neous case, however, the scalar field enters a regime of ampli- 
fied oscillations inside the overdensity after the turn around 
(as the (j)c term in Eq. 1461 becomes anti-frictional) , hence 
its oscillations are progressively larger. Consequently, the 
oscillations in 5c are more probable of becoming visible. 



4 SUMMARY AND CONCLUSIONS 

During the highly non-linear regime of cosmological pertur- 
bations, when dark matter halos are formed, scalar fields 
may become inhomogeneous. Hence, evolving differently in- 
side a collapsed overdensity and in the background Universe. 
The reason is simple: Different local spacetime geometries 
lead to different local equations of motion for the field. In 
the case of dark energy candidates, such as quintessence, 
this may result on dark energy having a different equation 
of motion in the background Universe and inside clusters. 
As a consequence, the dynamics of non-linear structure for- 
mation in the universe may show a distinct signature asso- 
ciated to the nature of dark energy and a particular model 
l|Mota fc van de Bruc3l2004^ . 

The behaviour of dark energy during the non-linear 
regime of structure formation can be divided into two ma- 
jor groups: the inhomogeneous models, which are the ones 
where the field inside the overdensity behaves differently 
form the background Universe and the homogeneous ones, 
where dark energy is just a smooth fluid which bath the 
whole Universe. 

In this article we have investigated the effects of inho- 
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mogeneities in dark energy on the time of collapse of a dark 
matter overdensity, and how this is reflected on the linearly 
extrapolated density threshold above which structures will 
end up collapsing, i.e. 5c{z) — Sl{z = Zcoi)- We have studied 
both coupled and minimally coupled scalar field candidates 
to dark energy as well as fluids with a constant equation 
of state. We have compared all these models and we found 
distinct features among them when taking into account the 
possible inhomogeneities in the dark energy fluid during the 
non- linear regime. 

In general Sc{z) depends on the particular model one 
considers. For constant equation of state candidates to dark 
energy there are minor differences between homogeneous 
and inhomogeneous models. The exception to that are the 
phantom energy candidates, with < —I. These present 
a quite distinct behaviour when one takes into account the 
inhomogeneities in the dark energy fluid. However, it is on 
scalar fields candidates where inhomogeneities in the dark 
energy fluid may result in more concrete features. The later 
will depend on its dynamical equation of state, which can 
be scale dependent, and on the contribution of the field to 
the energy budget of the Universe both at low and high 
redshifts. Coupled quintessence models behave similarly to 
the minimally coupled field candidates, but naturally show 
a strong dependence on the amount of dark matter coupled 
to the scalar field and the strength of the coupling. 

The reason for such a difference between inhomogeneous 
and homogeneous models can be understood looking at the 
Poisson equation 1251 which gives the evolution of the grav- 
itational potential inside the clusters. In opposition to the 
homogeneous models, one cannot neglect the importance of 
dark energy inside the matter overdensities. This leads to 
an extra source term in the usual Poisson equation, which 
affects the overall behaviour of the gravitational potential. 
This effect is specially important when the collapsing region 
enters into the non-linear regime. In particular, the time 
of collapse of the overdensity depends on the dark energy 
model one considers and on the equation of state inside it. 
These effects are reflected into linearly extrapolated density 
threshold above which structures will end up collapsing. Sc. 

The results of this work may have important con- 
sequences on cosmological studies of structure formation. 
In particular one can study its implications for the num- 
ber density of dark matter halos, their density profiles, 
galaxy number counting and dark matte r halo concentra- 
tions llNunes. da Silva fc AghanimI 120051: iManera fc Motal 
I2OO4I) . These investigations are imperative for cosmological 
studies that rely on these ingredients to measure dark en- 
ergy. Examples of these studies, for the case of homogeneous 
dark energy model s, include semi-analytical studies of strong 
lensing statistics iBartelmann et all |2003|: [Lopes fc Millei 
|20^) and weak lensing number counts jBMtehnMU^^T 
l2003f) . And for inhomogeneous mo dels include a study 
of da rk matter halo concentrations llLopes. Mota fc Milled 
12004) . 

Many other astrophysical phenomena may refiect the 
inhomogeneities in dark energy fluid during the non-linear 
regime. For instance, the merger of dark matter halos de- 
pend on the gravitational potential between them. Effects on 
the evolution and magnitude of the gravitational potential, 
may also affect peculiar velocities of galaxies inside clusters. 
One should also point out here, that in low-density regions 



of the universe, such as voids, inhomogeneous dark energy 
models, may play an important role as well. In those regions 
the local spacetime geometry may be different from the usual 
Friedman-Robertson- Walker metric one usually assumes for 
the background Universe. It is natural to think that dark en- 
ergy inside voids may behave differently from its background 
behaviour. This may again affect the dynamics of formation 
of voids in the universe. All these effects could in principle be 
detected using g alaxy redshift surveys such as the Sloan Dig - 
ital Sky Survey jAbazaiian et al.ll2003l: IZehavi et al]|2002l) . 
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